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Abstract 

Background: Smydl, the founding member of the Smyd family including Smyd-1, 2, 3, 4 and 5, is a SET and MYND domain 
containing protein that plays a key role in myofibril assembly in skeletal and cardiac muscles. Bioinformatic analysis revealed 
that zebrafish genome contains two highly related smydl genes, smydla and smydlb. Although Smydlb function is well 
characterized in skeletal and cardiac muscles, the function of Smydla is, however, unknown. 

Methodology/Principal Findings:lo investigate the function of Smydla in muscle development, we isolated smydla from 
zebrafish, and characterized its expression and function during muscle development via gene knockdown and transgenic 
expression approaches. The results showed that smydla was strongly expressed in skeletal muscles of zebrafish embryos. 
Functional analysis revealed that knockdown of smydla alone had no significant effect on myofibril assembly in zebrafish 
skeletal muscles. However, knockdown of smydla and smydlb together resulted in a complete disruption of myofibril 
organization in skeletal muscles, a phenotype stronger than knockdown of smydla or smydlb alone. Moreover, ectopic 
expression of zebrafish smydla or mouse Smydl transgene could rescue the myofibril defects from the smydlb knockdown 
in zebrafish embryos. 

Conclusion/Significance: Collectively, these data indicate that Smydla and Smydlb share similar biological activity in 
myofibril assembly in zebrafish embryos. However, Smydlb appears to play a major role in this process. 
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Introduction 

Members of the Smyd family are newly identified proteins that 
have been implicated in diverse biological functions in embryonic 
development and cancer [1]. Currently, five smyd genes (Smydl, 
— 2, —3, —4, and —5) have been identified in vertebrates based on 
the presence of both SET and MYND domains in their protein 
sequences [1]. Smydl, also known as skm-Bop, represents the first 
identified member of the Smyd family [2,3]. Smydl is specifically 
expressed in skeletal and cardiac muscles and plays a key role in 
muscle development and embryonic survival in mice and zebrafish 
[4—6]. Targeted disruption of the smydl gene resulted in defective 
cardiac morphogenesis and early embryonic lethality of mouse 
embryos [4] . Knockdown or mutation of smydlb gene in zebrafish 
led to disruption of myofibril organization in skeletal and cardiac 
muscles in zebrafish embryos [5,6]. 

The smydl gene is a direct downstream gene target of myogenic 
regulatory factors MyoD, Myogenin and Mef2 that control the 
muscle specific expression of smydl in skeletal muscles during 
embryogenesis and in adult muscle tissues [7-10]. A recent report 
showed that smydl expression is also regulated by serum response 
factor (SRF) through direct binding to the promoter region of 



smydl [11]. In addition, smydl gene expression can be repressed by 
Hepatoma-derived growth factor through interaction with a 
transcriptional co-repressor C-terminal binding protein (CtBP) 
[12]. Consistent with the idea of being a downstream factor of 
MyoD and Mef2, loss of Smyd 1 function had no effect on myoD 
and myogenin gene expression and myoblast specification [5]. 
However, loss of Smydl function resulted in defective sarcomere 
organization in myofibers of skeletal and cardiac muscles, 
suggesting that Smyd 1 is required in the late stage of muscle cell 
differentiation and myofiber maturation [5,6]. 

At present, little is known about the mechanism by which 
Smydl functions in myofibrillogenesis. In vitro studies have shown 
that Smydl has a histone methyltransferase (HMTase) activity 
[5,13], and could function as a transcriptional repressor in a 
histone deacetylase (HDAC)-dependent manor [4,14]. However, 
Just and colleagues reported recendy that the Smydl mutant 
lacking the HMTase activity was biologically active in myofibril 
assembly [6], arguing against Smydl being a HMTase activity 
dependent transcriptional regulator. Interestingly, Just and 
colleagues showed that GST-tagged Smydl was capable of pulling 
down skeletal muscle-specific myosin heavy chain [6]. Consistent 
with a potential role of Smyd 1 outside of the nucleus, a nuclear to 
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cytoplasmic translocation was observed during myoblast differen- 
tiation into myotubes [15], and Smydl is localized on the M-lines 
of sarcomeres although the biological significance of the sarco- 
meric localization is not clear [6,16]. 

Recent studies demonstrated that zebrafish genome contains 
two highly related smydl genes, smydla and smydl b [17]. Most of 
the previous studies were focused on smydlb, very little is known 
about smydla. It is not clear whether Smydla plays a similar role as 
Smydlb in myofibrillogenesis. To investigate the function of 
Smydla in myofibril assembly, we isolated smydla from zebrafish, 
and characterized its expression and function during muscle 
development. The results showed that smydla was specifically 
expressed in skeletal muscles of zebrafish embryos, smydla 
expression came several hours later than smydlb during myogenesis 
in zebrafish embryos. Functional analysis revealed that knockdown 
of smydla alone had little effect on myofibril assembly in zebrafish 
skeletal muscles. However, knockdown of smydla and smydlb 
together resulted in a stronger phenotype in myofibril disorgani- 
zation. Moreover, the myofibril defects from smydlb knockdown 
could be rescued by an ectopic expression of the zebrafish smydla 
or mouse Smydl transgene. Together, these data indicate that 
Smydla and Smydlb share similar biological activity in myofibril 
assembly although the function of Smydlb appears to be more 
critical. 

Results 

1. Characterization of Smydla in Zebrafish 

Sequence analysis revealed that zebrafish genome contains two 
highly related smydl genes (smydla and smydlb) with similar gene 
structures [17]. smydla and smydlb are believed to be generated by 
gene duplication. The zebrafish smydla is located on chromosome 
5, whereas smydlb is located on chromosome 8 (Figure SI). 
Sequence analysis revealed a strong synteny arrangement of 
zebrafish smydla gene and human smydl gene with the thronine 
synthase like 2 gene (THNSL2) and fatty acid binding proteinl gene 
(FABP1) in zebrafish and mouse genome (Figure SI). A similar 
synteny arrangement was found with the zebrafish smydlb and fatty 
acid binding protein lb like gene fabplb). However, no synteny 
arrangement was found with the smydlb gene and thronine synthase 
like 2 gene on chromosome 8 in zebrafish. Sequence alignment 
revealed that zebrafish Smyd 1 a contains the highly conserved SET 
and MYND domains involved in protein methylation and protein- 
protein interactions, respectively (Figure SI). Both Smydla and 
Smydlb share high sequence identity with Smydl from other 
vertebrate species (Figure S2), although only smydl gene has been 
identified in mice and humans. 

2. Temporal and Spatial Expression of Smydla in 
Zebrafish Embryos 

It has been reported previously that both smydla and smydlb are 
specifically expressed in developing muscles of zebrafish embryos 
[5,9,17]. To compare their patterns of expression during 
embryonic development, we carried out a temporal and spatial 
expression analysis by RT-PCR and whole mount in situ 
hybridization in zebrafish embryos. Data from the expression 
analysis showed that smydla and smydlb had different pattern of 
temporal expression, smydla was expressed several hours later than 
smydlb in zebrafish embryos, smydlb expression was first detected 
at 6 hours post-fertilization (hp!) with a strong expression starting 
around 9 hpf [5]. In contrast, smydla expression could not be 
detected until 19 hpf (Fig. 1A). 

Previous studies have demonstrated that smydlb encodes two 
muscle-specific isoforms of mRNA transcripts, designated as 



smydlb_Xw\ and smydlb_tw2 [5,16]. smydlb_tvl and smydlbjcv2 
are generated by alternative splicing of the 39 bp exon 5. 
smydlb_Xw\ transcript contains the 39 bp sequence from exon 5, 
whereas smydl b_tv2 transcript lacks this 39 bp sequence and thus 
is 13 aa shorter than smydlb_tvl [5,16]. Sequence analysis 
revealed that the zebrafish smydla gene also contained the small 
exon 5 of 39 bp. To test whether exon 5 could be alternatively 
spliced in smydla, we performed RT-PCR analysis using cDNA 
from zebrafish larvae of different stages. In contrast to smydlb 
where two different isoforms were amplified (Fig. IB), only the 
longer isoform of smydla was found to be expressed (Fig. 1C), 
suggesting that exon 5 in smydla was not alternatively spliced in 
zebrafish embryos. 

To determine whether smydla expression is restricted to muscle 
cells, we analyzed its spatial pattern of expression by whole mount 
in situ hybridization. The data showed that similar to smydlb, 
smydla is specifically expressed in skeletal muscles of zebrafish 
embryos (Fig. 1D-I). However, unlike smydlb, little or no 
expression was detected for smydla transcripts in cardiac muscles 
(Fig. ID, G). 

3. Knockdown of smydla Expression in Zebrafish 
Embryos had Little Effect on Muscle Development 

To determine whether smydla plays a role in muscle develop- 
ment, we knocked down smydla expression in zebrafish embryos 
using two smyd la-specific antisense morpholino oligos, E8I8-MO 
and E9I9-MO. The E8I8-MO and E9I9-MO were specifically 
targeted to the sequences at the exon-8/intron-8, or exon-9/ 
intron-9 junctions, respectively (Fig. 2A). RT-PCR analysis 
confirmed that both E8I8-MO and E9I9-MO could knock down 
the normal splicing of smydla transcripts and resulted in the 
production of defectively spliced smydla mRNA (Fig. 2B). 
However, the knockdown was not complete with an efficiency of 
approximately 50% (Fig. 2B). Sequence analyses revealed that 
E8I8-MO caused a defective splicing at the exon 8 and intron 8 
junction, resulting in a 14 bp deletion at the end of exon 8 
(Fig. 2A). This 14 deletion caused a reading frame shift, leading to 
the production of a mutant protein without the 105 aa C-terminal 
sequence. Similarly, the E9I9-MO caused a defective splicing at 
the end of exon 9 with a 58 bp deletion (Fig. 2A), which resulted 
in the production of a mutant protein without the last 63 aa C- 
terminal sequence. 

The smydla knockdown embryos were examined morphologi- 
cally for several days following the MO injection. Unlike the 
smydlb knockdown embryos that had no skeletal and cardiac 
muscle contraction, the smydla knockdown embryos had normal 
skeletal and cardiac muscle contraction when observed at 24, 48 
and 72 hpf. The E8I8-MO injected embryos appeared morpho- 
logically normal (Fig. 2E, F) compared with the control (Fig. 2C, 
D). Embryos injected with the smydla E9I9-MO, however, showed 
small head and developmental delay, suggesting some off-target 
toxic effect. The off-target toxic effect could be alleviated by co- 
injection with the p53-MO (Fig. 2G, H). To confirm that 
knockdown of Smydla did not affect early muscle development, 
we analyzed the specification of slow and fast muscle precursors 
and their subsequent migration and differentiation. Compared 
with control (Fig. 21), MyoD expression appeared normal in slow 
and fast muscle precursors of smydla knockdown embryos (Fig. 21, 
J). In addition, fiber-type-specific expression of slow or fast myosin 
heavy chain (MHC) genes also appeared normal (Fig. 2M, N). 
Slow fibers were clearly localized at the superficial layer of the 
myotome in smydla knockdown embryos (Fig. 2M), suggesting that 
migration of slow muscle cells was not affected by smydla 
knockdown. 
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Figure 1. The temporal and spatial pattern of smydla expression in zebrafish embryos. A. RT-PCR results show temporal expression of 
smydla in zebrafish embryos from fertilization to day 5. Elongation factor 1-alpha (ef-la) was used as control. B. RT-PCR analysis shows the alternative 
splicing of smydlb exon 5 generating two isoforms of smydlb, smyd1b_M and smyd1b_tv2. C. RT-PCR analysis shows the lack of alternative splicing 
of exon 5 in smydla in zebrafish embryos. D-l. Whole mount in situ hybridization shows the spatial pattern of smydla mRNA expression using a dig- 
labeled antisense probe, smydla expression was detected in skeletal muscles of zebrafish embryos at 24 (D, E, F) and 48 (G, H. I) hpf. D, G represent 
the side view; E, H repre4sent the dorsal view; F, I represent the cross sections. 
doi:1 0.1 371 /journal.pone.0086808.g001 
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Figure 2. Knockdown of smydla splicing by E8I8-MO and E9I9-MO. A. Location of the E8I8-MO and E9I9-MO splicing blockers at the junction 
of exon8/intron8 and exon 9/intron 9, respectively. Defective splicing using a cryptic splicing site in exon 8 in E8I8-MO injected embryos resulted in a 
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deletion of 14. Defective splicing using a cryptic splicing site in exon 9 in E9I9-MO injected embryos resulted in a deletion of 58 bp. B. RT-PCR 
showing the defective splicing induced by the E8I8-MO or E9I9-MO or both MOs. Compared with the PCR results from the wild type (wt) control 
embryos where a single band was generated, two bands were observed in the E8I8-MO injected embryos. Sequence analysis revealed that the 
smaller band contained a 14 bp deletion at the end of exon 8. The upper band represented the heteroduplex formed by the normal and defectively 
spliced products. Similarly, three bands were detected in E919-MO injected embryos. Sequence analysis revealed that the smaller band contained a 
58 bp deletion at the end of exon 9. The middle band was the normal spliced product, whereas the upper band represented the heteroduplex 
formed by the normal and defectively spliced products. Two major bands were detected in E8I8-MO and E9I9-MO co-injected embryos. Sequence 
analysis revealed that the smaller band resulted from defective splicing. C-H. Morphology of smydla knockdown embryos at 48 hpf. Zebrafish 
embryos injected with control MO (C, D), E8I8-MO (E, F) or E9I9-MO (G, H) at 48 hpf. I-N. Knockdown of smydla expression had no effect on myoblast 
specification and early differentiation of slow and fast myofibers. In situ hybridization showing normal MyoD expression in control-MO (I) or smydla 
MO (L) injected embryos at 14 hpf. Cross views of slow (J, M) or fast (K, N) MHC expression in control (J, K) or smydla MO (M, N) injected embryos at 
24 hpf. 

doi:1 0.1 371 /journal.pone.0086808.g002 



4. Knockdown of smydla and smydlb together Resulted 
in a Stronger Muscle Phenotype in Zebrafish Embryos 

To analyze the effects of smydla knockdown on sarcomere 
assembly in skeletal muscles, we first examined the organization of 
thick filaments in slow myofibers of smydla knockdown zebrafish 
embryos. The results showed that knockdown of smydla alone had 
no detectable effect on the sarcomere organization of the myosin 
thick filaments (Fig. 3D-F) in slow fibers of zebrafish embryos at 
24, 48 and 72 hpf. The lack of muscle defect from smydla 
knockdown could be due to the inefficient knockdown. To 
determine whether increasing the efficiency of smydla knockdown 
could affect muscle development, we co-injected the two smydla 
splicing MOs into zebrafish embryos. Compared with the single 
MO injection, the results showed that co-injection of the two MOs 
increased the efficiency of smydla knockdown (Fig. 2B). However, 
it did not result in more defects on myofibril assembly in slow 
muscles (Fig. 3G-I). 

To test the possibility that the lack of muscle defect from smydla 
knockdown was due to the redundant function from smydlb, we 
carried out the double knockdown of smydla and smydlb in 
zebrafish embryos. The data showed that knockdown of smydla 
and smydlb together resulted in a stronger muscle phenotype 
compared with smydla or smydlb knockdown alone. As shown in 
Fig. 3J, knockdown of smydlb alone resulted in disruption of 
myosin thick filament organization at 28 hpf. However, the 
defective thick filament organization was partially recovered in 
smydlb knockdown zebrafish embryos at 48 and 72 hpf (Fig. 3K, 
L). Double knockdown studies showed that the partial recovery 
from smydlb knockdown was diminished when smydlb was knocked 
down together with smydla (Fig. 3N,0), suggesting that Smydla 
may have a partial function redundancy with Smydlb in myofibril 
assembly. 

To determine whether this was also true for other sarcomeric 
structures, we analyzed the ot-actin thin filaments and Z-line 
organization in the single or double knockdown embryos. The 
results showed that the knockdown of smydla and smydlb together 
had a stronger effect on thin filament disruption at 72 h pf 
(Fig. 40) compared with the knockdown alone (Fig. 41, L). A 
similar finding was also observed with the Z-line structure (Fig. 5). 
Collectively, these data indicate that Smydla and Smydlb are 
required for myofibril organization and assembly in skeletal 
muscles of zebrafish embryos although Smydlb plays a dominant 
role compared with Smydla. 

In situ expression analysis showed that smydla was expressed in 
fast muscles (Fig. 1). To determine whether knockdown of smydla 
had any effect on fast muscles, sarcomere organization was 
characterized in fast muscles of smydla knockdown embryos. The 
results showed that knockdown of smydla alone had no effect on 
myofibril assembly in fast muscles (Fig. 6). 



5. Ectopic Expression of Zebrafish smydla or Mouse 
Smydl Transgene Could Rescue the Myofibril Defects in 
smydlb Knockdown Zebrafish Embryos 

It has been suggested that actinopterygian fish, such as 
zebrafish, had gone through one more round of whole genome 
duplication compared with mammals during evolution [18]. The 
two Smydl genes in zebrafish were likely generated by the gene 
duplication event in teleosts because mouse and human genomes 
contain only one Smydl gene. Sequence analysis revealed that 
zebrafish Smyd 1 a and Smyd 1 b share high sequence similarity with 
the mouse and human Smydl (Fig. S2). However, functional 
studies revealed that knockdown of smydla or smydlb resulted in 
strong differences in muscle phenotype. Knockdown of smydlb 
significantly disrupted the myofibril organization, whereas knock- 
down of smydla alone had little effect. The different severity of the 
knockdown phenotypes could be due to structural changes of 
smydla protein that made it less active compared with smydlb. 
Alternatively, the different knockdown phenotype could be caused 
by their distinct patterns of temporal expression because smydla 
was expressed several hours later than smydlb in zebrafish 
embryos. 

To clarify these questions, we tested whether ectopic expression 
of smydla early in zebrafish embryos could rescue the smydlb 
knockdown phenotypes. A rescue experiment was performed in 
the smydlb knockdown zebrafish embryos using the smydla 
transgene (Smydl -zfsmydl a my ") directed by the smydlb promoter. 
The Smydl -zfsmydl a myc transgene expressing a myc-tagged zebra- 
fish Smydla was co-injected with smydlb ATG-MO into zebrafish 
embryos. The smydlb ATG-MO could specifically knock down the 
expression of smydlb, however it had no effect on the expression of 
the zebrafish smydla™ 3 " transgene because the transgene did not 
contain the smydlb ATG-MO target sequence. Myofibril organi- 
zation and ectopic Smydla my< expression was analyzed in the co- 
injected zebrafish embryos by double staining with anti-MyHC 
(F59) and an anti-myc antibodies. The data showed a clear rescue 
of thick filament organization in the co-injected embryos (Fig. 7A). 
The rescue appeared in a mosaic pattern, consistent with the 
pattern of gene expression through DNA injection. Double 
staining revealed a perfect match between the rescued myofibers 
and the expression of the myc-tagged zebrafish Smydla (Fig. 7C, 
E). Collectively, these data indicate that Smydla could replace 
Smydlb function if expressed early in muscle cells, arguing that 
the different severity of smydla and smydlb knockdown phenotypes 
was likely due to distinct pattern of temporal expression. 

Sequence analysis revealed that Smydla and Smydlb share 
similar sequence identity with mouse and human Smydl (Fig. S2). 
To test whether Smydl function is conserved during evolution, we 
performed a rescue experiment in the smydlb knockdown zebrafish 
embryos using the mouse smydl transgene directed by the zebrafish 
smydlb promoter. The mouse Smydl transgene (smydl -mSmydl myc ) 
expressing a myc-tagged mouse Smydl was co-injected with the 
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Figure 3. The effect of smydla and/or smydlb knockdown on myosin thick filament organization in skeletal muscles. Zebrafish 
embryos injected with smydlb MO or smydla MO or both were fixed at 28, 48 and 72 hpf. Myosin thick filament organization was analyzed by 
immunostaining with the F59 antibody which recognizes the myosin heavy chain (MHC) in slow muscles, and followed by FTIC-labeled secondary 
antibody. The images represent side view of trunk muscles around segment 1 0. A-C. Lateral view of thick filament organization in slow muscle fibers 
of control-MO injected embryos at 28 (A), 48 (B) and 72 (C) hpf. D-F. Lateral view of thick filament organization in slow muscle fibers of smydla E8I8- 
MO injected embryos at 28 (D), 48 (E) and 72 (F) hpf. G-l. Lateral view of thick filament organization in slow muscle fibers of smydla E8I8-MO and 
E9I9-MO co-injected embryos at 28 (G), 48 (H) and 72 (I) hpf. J-L. Lateral view of thick filament organization in slow muscle fibers of smydlb ATG-MO 
injected embryos at 28 (J), 48 (K) and 72 (L) hpf. M-O. Lateral view of thick filament organization in slow muscle fibers of smydla E8I8-MO and 
smydlb ATG-MO co-injected embryos at 28 (M), 48 (N) and 72 (0) hpf. Scale bars: 20 urn in A-C. 
doi:1 0.1 371 /journal.pone.0086808.g003 



smydlb ATG-MO into zebrafish embryos. Double staining 
revealed that expression of the myc-tagged mouse Smydl could 
rescue the myofibril defects from the smydlb knockdown in a cell 
autonomous manner (Fig. 7B, D, F). Collectively, these data 
indicate that Smydl function in myofibril organization is likely 
conserved during evolution. 

6. The Subcellular Localization of Zebrafish Smydla and 
Mouse Smydl in Myofibers 

Previous studies by us and others have shown that Smydlb is 
localized on the M-line of sarcomeres in both skeletal and cardiac 
muscles [6,16]. This was especially evident for the longer isoform 
Smydl b_tvl [16]. It is not clear whether zebrafish Smydla and 
mouse Smydl could also localize on M-lines of sarcomere in 
skeletal muscles. To clarify this question, we examined the 
subcellular localization of myc-tagged zebrafish Smydla and 



mouse Smydl in skeletal muscles of zebrafish embryos. The results 
showed that, similar to Smydlb, mouse Smydl was localized on 
the sarcomere of skeletal muscles (Fig. 7D). However, unlike 
Smydlb, zebrafish Smydla showed no clear sarcomeric localiza- 
tion in myofibers (Fig. 7C). 

To enhance the sensitivity of the detection, we generated a 
transgene, pTol2-Smydla -EGFP, expressing a GFP-tagged zebra- 
fish Smydla fusion protein. The pTol2-Smydl a-EGFP transgene 
was injected into zebrafish embryos. Smydl a-EGFP expression 
and subcellular localization were characterized in the injected 
zebrafish embryos by direct observation. The results showed that 
Smydl b-EGFP had a clear sarcomeric localization (Fig. 8B, E). In 
contrast, Smydl a-EGFP was not localized on the sarcomeres 
(Fig. 8 A, D). A diffuse cytosolic localization was detected for 
Smydl a-EGFP in myofibers of zebarfish embryos (Fig. 8A, B), 
similar to EGFP control (Fig. 8C, F). Together, these data indicate 
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Figure 4. The effect of smydla or smydlb single or double knockdown on thin filament organization in skeletal muscles. Zebrafish 
embryos injected with smydlb MO or smydla MO or both were fixed at 28, 48 and 72 hpf. a-actin thin filament organization was analyzed by 
immunostaining with anti-a-actin antibody (Acl-20.4.2), and followed by FTIC-labeled secondary antibody. The images represent side view of trunk 
muscles around segment 10. A-C. Lateral view of thin filament organization in skeletal muscle fibers of control-MO injected embryos at 28 (A), 48 (B) 
and 72 (C) hpf. D-F. Lateral view of thin filament organization in skeletal muscle fibers of smydla E8I8-MO injected embryos at 28 (D), 48 (E) and 72 
(F) hpf. G-l. Lateral view of thin filament organization in slow muscle fibers of smydla E8I8-MO and E9I9-MO co-injected embryos at 28 (G), 48 (H) and 
72 (I) hpf. J-L. Lateral view of thin filament organization in skeletal muscle fibers of smydlb ATG-MO injected embryos at 28 (J), 48 (K) and 72 (L) hpf. 
M-O. Lateral view of thin filament organization in skeletal muscle fibers of smydla E8I8-MO and smydlb ATG-MO co-injected embryos at 28 (M), 48 
(N) and 72 (O) hpf. Scale bars: 20 urn in A-C. 
doi:1 0.1 371 /journal.pone.0086808.g004 



that Smydla and Smydlb have distinct subcellular localization, 
indicating that the sarcomeric localization of Smyd 1 b might not be 
directly linked with the Smyd 1 function in sarcomere organization 
because Smydla could rescue the myofibril defects from Smydlb 
knockdown although Smydla was not localized on the sarcomere. 
This is consistent with previous findings that either one of the 
Smydlb isoforms, Smydlb_tvl or Smydlb_tv2, from alternative 
splicing was able to rescue the smydlb knockdown phenotype 
although their proteins showed distinct subcellular localization 
[5,16].' 

Discussion 

In this study, we have characterized the muscle phenotype from 
knockdown of smydla or smydlb alone, or together in zebrafish 
embryos. We demonstrated that in contrast to smydlb which is 
absolutely required for myofibril organization in skeletal muscles, 



knockdown of smydla alone had very littie effect on myofibril 
organization. However, knockdown of smydla and smydlb together 
resulted in a stronger muscle phenotype compared with knock- 
down of smydla or smydlb alone. We further demonstrated that the 
muscle defects from smydlb knockdown could be rescued by the 
ectopic expression of the zebrafish smydla or mouse Smydl, 
suggesting that Smydl function is likely to be conserved during 
evolution. Finally, we showed that similar to Smydlb, mouse 
Smydl was localized on sarcomeres in skeletal myofibers. In 
contrast, zebrafish Smydla did not show any sarcomeric 
localization although it could rescue the Smydlb knockdown 
defects. 
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Figure 5. The effect of smydla or smydlb single or double knockdown on the Z-line organization in skeletal muscles. Zebrafish 
embryos injected with smydlb MO or smydla MO or both were fixed at 28, 48 and 72 hpf. Z-line organization was analyzed by immunostaining with 
anti-a-actinin antibody (EA-53), and followed by FTIC-labeled secondary antibody. The images represent side view of trunk muscles around segment 
10. A-C. Lateral view of Z-line organization in skeletal muscle fibers of control-MO injected embryos at 28 (A), 48 (B) and 72 (C) hpf. D-F. Lateral view 
of Z-line organization in skeletal muscle fibers of smydla E8I8-MO injected embryos at 28 (D), 48 (E) and 72 (F) hpf. G-l. Lateral view of Z-line 
organization in skeletal muscle fibers of smydla E8I8-MO and E9I9-MO co-injected embryos at 28 (G), 48 (H) and 72 (I) hpf. J-L. Lateral view of Z-line 
organization in skeletal muscle fibers of smydlb ATG-MO injected embryos at 28 (J), 48 (K) and 72 (L) hpf. Wl-O. Lateral view of Z-line organization in 
skeletal muscle fibers of smydla E8I8-MO and smydlb ATG-MO co-injected embryos at 28 (M), 48 (N) and 72 (O) hpf. Scale bars: 20 urn in A-C. 
doi:1 0.1 371 /journal.pone.0086808.g005 



The Unequal Role of Smydla and Smydlb in Myofibril 
Assembly 

Gene duplication plays a vital role in evolution [19]. Sequence 
analysis revealed that zebrafish genome contains two smjdl genes, 
smydla and smydlb that were likely generated by the whole genome 
duplication in actinopterygian fish during evolution [17,18,20]. 
smydla and smydlb are paralogous genes that are dispersed on 
different chromosomes in zebrafish. Sequence analysis revealed 
that zebrafsh smydla and mouse Smydl showed a highly conserved 
synteny arrangement with the surrounding THNSL2 and FABP1 
genes. The synteny arrangement with fabplwas maintained for 
smydlb in zebrafish. However, this synteny was not found between 
the zebrafish smydlb and thnsl2 genes, suggesting that a genomic 
rearrangement or deletion might occur in this region after the 
whole genome duplication. 

It has been suggested that after whole genome duplication, the 
duplicated genes could evolve in several directions without the 
selective pressure. Some of the duplicated genes could gain new 



functions, while many others turned into nonfunctional pseudo- 
genes [19,20]. We showed in this study that both smydla and 
smydlb are functional genes after the duplication. However, they 
do not play an equal role in myofibril organization. Smydlb 
appears to be more critical than Smydla because knockdown of 
smydlb gave a stronger muscle phenotype compared with smydla 
knockdown. Our data indicate that functional difference between 
Smydla and Smydlb could be caused by their different patterns of 
temporal expression, different activity due to changes of protein 
structure, or different levels of gene expression. We showed that 
early ectopic expression of smydla by using the smydlb promoter 
could rescue the myofibril defects from smydlb knockdown, 
suggesting that Smydla could function as Smydlb if expressed 
in a similar temporal and spatial pattern in zebrafish embryos. At 
present, it is not clear whether Smydla gained new functions 
during evolution. Knockdown of smydla alone failed to reveal any 
significant phenotype in zebrafish embryos, unless it was knocked 
down together with smydlb. Together, these studies suggest that 
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Figure 6. The effect of smydla knockdown on the sarcomere organization in fast muscles. Zebrafish embryos co-injected with smydla 
E8I8 and E9I9 MOs were fixed at 28, 48 and 72 hpf. The sarcomere organization of thin filaments and Z-lines was analyzed by immunostaining with 
anti-a-actin and anti-a-actinin antibody, respectively, and followed by FTIC-labeled secondary antibody. The images represent side view of trunk 
muscles around segment 10. A-C. Lateral view of thin filament organization in skeletal muscle fibers of control-MO injected embryos at 28 (A), 48 (B) 
and 72 (C) hpf. D-F. Lateral view of thin filament organization in skeletal muscle fibers of smydla MO injected embryos at 28 (D), 48 (E) and 72 (F) hpf. 
G-l. Lateral view of Z-line organization in skeletal muscle fibers of control-MO injected embryos at 28 (G), 48 (H) and 72 (I) hpf. J-L. Lateral view of 
the Z-line organization in skeletal muscle fibers of smydla MO injected embryos at 28 (J), 48 (K) and 72 (L) hpf. Scale bars: 20 u,m in D-F and J-L. 
doi:1 0.1 371 /journal.pone.0086808.g006 



Smyd 1 a is involved in myofibril assembly, but with a minor role in 
the process. 

The SmycM Function in Myofibril Organization is 
Conserved during Evolution 

Our studies demonstrated that Smydl function in myofibril 
organization is likely to be conserved during evolution because 
ectopic expression of mouse Smydl could rescue the muscle defects 
from smjdlb knockdown in zebrafish embryos. Unlike zebrafish, 
there is only one Smydl gene in mouse and human genomes. 
Consistent with the idea that zebrafish smjdlb is the ortholog of 
Smydl in mice, it has been reported that knockout of Smydl in mice 
resulted in early embryonic lethality from cardiomyogenesis 
defects [4]. A cardiac muscle defect was also observed in Smyd lb 
knockdown or mutant zebrafish embryos [5,6] . In addition to their 
functional similarity, the zebrafish Smyd lb and mouse Smydl 
share other common features in gene expression and protein 
localization. First, both zebrafish smydl b and mouse smydl express 
two alternatively spliced mRNA transcripts that differ by 39 bp 
encoded by exon 5. Second, both mouse Smydl and zebrafsh 
Smyd 1 b are localized on the M-lines of sarcomeres in skeletal and 
cardiac muscles [6], which is in contrast to the zebrafish Smydla 
that showed no sarcomeric localization. 



The Sarcomeric Localization of Smydl 

The biological significance of the sarcomeric localization is not 
clear. It has been reported that Smyd lb could bind to myosin and 
showed a transient localization with thick filaments during 
myofibrillogenesis in zebrafish embryos [6]. However, after the 
completion of sarcomere organization, the Smyd 1 b is translocated 
to the M-line. Interesting, Etard and colleagues have shown that 
myosin chaperones Unc45b and Hsp90ot could shuttie between 
the A band and the Z line in response to stress or damage to the 
myofiber [21]. It has been suggested that the sarcomeric 
localization of myosin chaperones and their translocation within 
different parts of the muscle cells could be involved in the response 
of muscle cells to mount efficient physiological responses to muscle 
stress, load requirements, and/or stretch [21]. It remains to be 
determined whether the sarcomeric localization of Smyd lb is 
involved in sarcomere remodeling and response to muscle stress. 

We have shown previously that the sarcomeric localization of 
Smydlb requires Phe223 and Serine 225 [16]. Substitution of 
Phe223 or Ser225 with alanine significantly diminished the 
sarcomeric localization of Smydlb [16]. Interestingly, the two 
highly conserved residues Phe223 and Ser225 are present in 
zebrafish Smydla, arguing that in addition to the Phe223 and 
Se225 residues, other residues or motifs are likely to be involved in 
the sarcomeric localization. Sequence analysis revealed that 
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Figure 7. The myofibril defects from smydlb knockdown could be rescued by ectopic expression of zebrafish smydla or mouse 
Smydl. The smydlb ATG-MO was co-injected with smydl b:zfsmyd1 a myc or smydl b:mSmydl myc transgene into zebrafish embryos at 1 or 2 cells 
stages. Myosin thick filament organization and transgene expression were analyzed by double immunostaining with anti-MHC (F59, green) and anti- 
myc (9E10, red) antibodies. A, C, E. Double immunostaining with anti-MHC (A) and anti-myc (C) antibodies shows the normal thick filament 
organization in myofibers expressing the myc-tagged zebrafish smydla transgene at 24 hpf. E, merged picture of A and C. B, D, F. Double 
immunostaining with anti-MHC (B) and anti-myc (D) antibodies shows the normal thick filament organization in myofibers expressing the myc- 
tagged mouse Smydl transgene at 24 hpf. F, merged picture of B with D. D and F showed the sarcomeric localization of myc-tagged mouse Smydl 
(red). The myc-tagged mouse Smydl was localized in the middle of the MHC thick filament (F), a region normally occupied by the M-line. Scale bars: 
20 urn in E; 14 u,m in F. 
doi:1 0.1 371 /journal.pone.0086808.g007 



Smydla differs from Smydlb significantly at the C-terminal 
sequence (Figure SI). Compared with Smydla, Smydlb has seven 
extra amino acids at the C-terminus including three highly 
conserved residues (xxLFxxK) that are present in all Smydl 
proteins that have been identified in vertebrates. It remains to be 
determined whether these residues may contribute to the 
sarcomeric localization of Smydlb. 

Materials and Methods 

Ethics Statement 

This study was carried out in strict accordance with the 
recommendations in the Guide for the Care and Use of 
Laboratory Animals of the National Institutes of Health. The 
protocol was approved by the Institutional Animal Care and Use 
Committee of the University of Maryland (Permit Num- 
ber:06 10009). 

Fish Maintenance and Use 

Zebrafish were maintained at 28.5 C in 10 gallon aquarium 
supplied with fresh water and air at a photoperiod of 14 hours of 



light and 10 hours of dark at the Zebrafish Facility at Aquaculture 
Research Center (ARC) in Columbus Center (Baltimore). Pairs of 
adult male and female zebrafish were put in 1 liter tanks for 
natural spawning after lights were turned on in the facility. 
Embryos were collected and raised at 28.5 C. To ease pain and 
facilitate animal handling, fish embryos over one day old were 
anaesthetized in 0.6 mM Tricaine that has been buffered to 
neutral pH around 7. The anaesthetized embryos were sacrificed 
quickly by keeping on ice and used directly for RNA extraction or 
fixation for immunostaining. All procedures were carried out in 
compliance with the guidelines stipulated by the Institutional 
Animal Care and Use Committee of the University of Maryland. 

Isolation of Zebrafish smydla cDNA 

The full-length coding sequence of zebrafish smydla 
(NM_205540) was cloned by RT-PCR using total RNA extracted 
from zebrafish embryos at 24-72 hpf. The PCR was carried with 
the Advantage PCR kit using the following primers: zismydla-Pl 
(5' -AGCATGACCGTGGAGAAGACGGAC-3 ') and zismydla- 
P2 (5' -GTGTTCATGCTTTGATCTGCACTT-3') derived 
from DNA sequences near the start codon and stop codon, 
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Figure 8. Characterization of the Smydla subcellular localization using the Smydla-EGFP fusion protein. DNA constructs expressing 
the zebrafish Smydla-EGFP or Smydlb-EGFP fusion proteins or EGFP control were injected into zebrafish embryos. Their expression and localization 
were determined in myofibers of the injected zebrafish embryos at 27 (A, B, C) and 48 (D, E, F) hpf. A and D represent Smydlb-EGFP; B and E 
represent Smydla-EGFP; C and F represent EGFP control. Scale bars: 6 urn in B; 10 um in E. 
doi:1 0.1 371 /journal.pone.0086808.g008 



respectively. The PCR product was cloned into pGEM-T vector. 
The resulting plasmid was named pGEM-T-zfrmyrfia. 

Analysis of smydla and Smdlb Expression in Zebrafish 
Embryos by RT-PCR 

To analyze the temporal pattern of smydla expression during 
development, total RNA was extracted from zebrafish embryos at 
20min, 3h, 6h, 9h, 12h, 16h, Id, 2d, 3d, 4d and 5d post 
fertilization, smydla expression was determined by RT-PCR using 
zismydla-Yl and zfsmydla-P2 primers. To characterize the 
expression of tvl and tv2 isoforms from alternative splicing of 
exon 5, RT-PCR was carried out using two new sets of primers 
smydla-tw-Pl/72 or smydlb-Xv-Pl/Y2, respectively. PCR using 
these primers produced small PCR fragments covering the 
alternatively spliced exon 5. 

smydla-w-P 1 :5 '-GTCTGAATCGGTCGGC ATCGGC-3 ' . 

jmy;fl-tv-P2:5'-CATAACTGACGGTCAGCTCCTGT-3'. 

smydl A-tv-P 1:5' -TCATGGTGAGCGATC AGC GC GGC-3 ' . 



imy/6-tv-P2:5'-ACGTTCAGATAATCCACATACGC-3'. 

Synthesis of Morpholino Antisense Oligos 

Morpholino antisense oligos were synthesized by Gene Tools. 
The smydlb ATG-MO was based on the target sequence 
surrounding the ATG start codon. The two smydla splicing 
blockers (E8I8-MO, E9I9-MO) were based on the sequence at the 
exon-8/intron-8 or exon-9/intron-9 junction, respectively. The 
standard control MO from Gene Tools was used as control. 

smydlb ATG-MO: 5' - ACTTCCACAAACTCCATTCTG- 
GATC-3'. 

smydla E8I8-MO: 5' - ATATCGCAACACTCACATG- 
TATCCA-3'. 

smydla E9I9-MO: 5' 
CAGGTCTCTGCTGA-3 ' . 



GGTTGTACCTC- 
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Morpholino and DNA Microinjection in Zebrafish 
Embryos 

Morpholino antisense oligos were dissolved in lx Danieau 
buffer [22] to a final concentration of 0.5 mM. Approximately 1— 
2 nl (5-10 ng) of MO was injected into each embryo at 1 or 2 cell 
stages. For co-injection with transgenes, equal volumes of MO 
(1 mM) and DNA construct (100 Hg/ml) was mixed and 1-2 nl 
was microinjected into each embryo at 1 or 2 cell stages. All the 
microinjection experiments and subsequent analyses were carried 
out at least three times with 100-150 embryos per sample in each 
experiment. 

Analysis of Smydla mRNA Splicing in MO Injected 
Embryos by RT-PCR 

To determine the efficacy of E8I8-MO and E9I9-MO splicing 
blocker on smydla transcript splicing, total RNA was extracted 
from E8I8-MO or E9I9-MO single injection or E8I8-MO and 
E9I9-MO co-injected embryos at 48 hpf. smydla transcripts were 
amplified using smydla-FLl-F\ and smydla-YI primers or smydla- 
E7/8-F1 and smydla-P2 primers. The smydl a-FH-YX primer (5' - 
CTTCAGGCTCTGGTGAAGATTGAA-3') was derived from 
the exon7. The smydl a-E7/8-F\ primer (5' -ACTTCCAT- 
GAGGTGATCAGGATCTG-3 ') was derived from the junction 
of the exon7 and exon 8 sequences, while the smydl a-F2 primer (5' 
-TCATGCTTTGATCTGCACTTG- 3') was from the antisense 
sequence near the stop codon. The PCR products were cloned 
into pGEM-T Easy vector for sequencing analysis. 

Construction of Zebrafish smydl b-zfsmyd1a myc 
Transgene 

To generate the smydl b-zfsmydla Tnyc transgene for muscle 
specific expression of myc-tagged zebrafish Smydla, the coding 
sequence of the zebrafish smydla was amplified with pfu DNA 
polymerase using the pGFM-T-zismydla plasmid as a template. 
The PCR was carried out using zismydla-FX and zfsmydla-'i' -my c 
primers. The zfsmydla-Pl primer (5 ' - AGC ATG AC C GTGG A- 
GAAGACGGAC- 3') contained the ATG start codon. The 
zfsmydl a-3'-myc primer (5'-CTAATTCAGGTCCTCTTCAGA- 
GATGAGCTTCTGCTCTGCTTTGATCTGCACTTGTTC- 
3') contained a myc-tag coding sequence and a stop codon. The 
PCR product containing the full-length coding sequence of a myc- 
tagged zebrafish smydla was phosphorylated with T4 kinase and 
directly cloned in the Smal sites of the muscle-specific smydl b-g/p 
vector [9]. The resulting plasmid was named .smydl b-zismydla Tayc . 

Construction of smydl b-mSmyd1 myc construct 

Mouse smydla cDNA was cloned by RT-PCR using total RNAs 
from adult mouse skeletal muscles. The PCR was carried with the 
Advantage PCR kit using the mouse mSmydl-ATG (BH1) and 
mSmydl-3'-myc primers. The PCR product was cloned into 
pGEM-T vector. The resulting plasmid was named pGEM-T- 
msmydla™ 7 ''. 

The mSmydl-ATG (BH1) primer (5'-GGATCCATGGA- 
GAACGTGGAGGTCTTC- 3') contained the ATG start codon, 
while the mSmydl-3'-myc primer (5'CTAATT- 
CAGGTCCTCTTCAGAGAT- 

GAGCTTCTGCTCCTGCTTCTTATGGAACAGAG-3') con- 
tained the myc-tag coding sequence and a stop codon. 

To generate the smydl b-mSmydl myl ' transgene for muscle specific 
expression of myc-tagged mouse Smydl, the coding sequence of a 
myc-tagged mouse Smydl was amplified from the pGEM-T- 
mSmydl myc plasmid by pfu DNA polymerase. The PCR was 
carried using the following primers: mouse Smydl-ATG (BH1) 



and mouse Smydl-3'-myc. The PCR product containing the 
coding sequence of mouse Smydl and the myc-tag at the C- 
terminus was phosphorylated with T4 kinase and directly cloned 
into the Smal sites of the muscle-specific smydl b vector [9]. The 
resulting plasmid was named smydl b-mSmydl myl . 

Construction of Tol2-zfsmyd1a-EGFP Construct 

The pTol2-OTp<</./<2-EGFP construct was generated by cloning 
the zebrafish smydla coding sequence in frame upstream of the 
EGFP coding sequence in the T2A200R150G vector [23]. Briefly, 
the zebrafish smydla coding sequence without the stop codon were 
generated by PCR using pfu DNA polymerase (Stratagene). The 
PCR was carried out using the smydlb-zfsmydla myc plasmid as a 
template together with the zismydla-F (5'- gctAGATCTatgaccgtg- 
gagaagacggaccc-3 ') and zismydla-H (5'- gctAGATCTgctttgatctg- 
cacttgttctcc-3') primers. A BgUI site was introduced at the 5' and 
3' ends of the smydla coding sequence via the respective PCR 
primers. The PCR products were digested with BgUI and then 
cloned into the compatible BamHI site of the T2A200R150G 
vector. The DNA sequence at the smydla and EGFP junction was 
confirmed by sequencing. 

Whole Mount in Situ Hybridization and Immunostaining 

Whole mount in situ hybridization was carried out using 
digoxigenin-labeled antisense probes as previously described [24— 
26]. The pGEM-T-zfswjx/ia plasmid was digested with Ncol and 
transcribed with Sp6 RNA polymerase to synthesize digoxigenin- 
labled antisense RNA probes. Immunostaining was carried out 
using whole mount zebrafish embryos (1-3 dpfj as previously 
described [5,27]. Briefly, zebarfish embryos were fixed in 4% 
paraformaldehyde (in PBS) for 1 hour at room temperature. The 
fixed embryos were washed for 15 minutes 3 times in PBST (0.5% 
Tween-20 in lx PBS). Embryos of 48 and 72 hpf were digested 
with 1 mg/ ml collagenase for 45 and 75 minutes, respectively. In 
situ hybridization for MyoD, slow and fast myosin heavy chains 
were carried out using the respective digoxigenin-labeled probes as 
descried previously [28]. 

Immunostaining was performed with the following primary 
antibodies: anti-a-actinin (clone EA-53, $A7811, Sigma), anti- 
MHC for slow muscles (F59, DSHB), anti-actin (Ac 1-20.4.2, 
Progen), and anti-myc tag (7 ID 10, Cell Signaling). Secondary 
antibodies were FTIC or TRIC-conjugates (Sigma). The images 
were acquired using an upright microscope (Zeiss, Oberkochen, 
Germany) equipped with a confocal image analyzer (BIO-RAD 
Radiance 2100 Imaging Systems, Hercules, CA). 

Supporting Information 

Figure SI Structure comparison of zebrafish Smydla 
and Smydlb. A. The synteny arrangement of zebrafish smydl b 
and mouse smydl with fatty acid binding protein lb (fabplb) and thronine 
synthase like 2 (thnsl2) genes. B. The diagrams of zebrafish smydla 
genomic and protein structures. C and D. Sequence comparison 
of functional MYND (C) and SET (D) domains in Smyd 1 a and 
Smydlb with other vertebrate Smydl and Smyd3 proteins. The 
identically conserved residues in the MYND and SET domains are 
indicated by the asterisk (*). The similar amino acid residues are 
indicated by = . ZF, zebrafish; Fu, fugu; Ch, chicken; Mo, mouse; 
Hu, human. 
(PDF) 

Figure S2 Sequence alignment of vertebrate Smydl 
proteins. Sequence comparison of zebrafish Smydla, Smydlb, 
chicken Smydl, mouse Smydl and human Smydl proteins. 
(PDF) 
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